Staggered magnetization, critical behavior and weak ferromagnetic properties of 

LaMnC>3 by muon spin rotation. 



M. Cestelli Guidi, G. Allodi, R. De Renzi, G. Guidi 
Dipartimento di Fisica e Istituto Nazionale di Fisica della Materia, 
Universita di Parma, F43100 Parma, Italy 



M. Hennion 

Laboratoire L. Brillouvn, CEA-CNRS, CE Saclay, F-91191 Gif sur Yvette 



France 



L. Pinsard 

Laboratoire de Chimie des Solides, Universite Paris-Sud, F-91405, France 



o 
o 

(N 



(N 

(N 
> 
oo 
in 

(N 

o 
o 

-I— > 



c 

o 
o 



X 



A. Amato 

Lab. for Muon-Spin Spectroscopy, Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland 

(Dated: February 1, 2008) 

We present a study of a microtwinned single crystal of LaMnOa by means of implanted muons. 
Two muon stopping sites are identified from the symmetry of the internal field in the ordered phase. 
The temperature dependence of these fields yields the behavior of the staggered magnetization 
from which a static critical exponent (/3 = 0.36(2)) is extracted and discussed. The muon spin- 
spin relaxation rate shows a critical slowing down (contrary to preliminary findings) with a critical 
exponent n = 0.7(1), witnessing the Ising nature of the dynamic fluctuations. The muon precession 
frequencies vs. applied magnetic field reveal the saturation of the weak ferromagnetic domain 
structure originated by the Dzialoshinski-Moriya antisymmetric exchange. 



I. INTRODUCTION 

LaMn03 is a layer antiferromagnetic perovskite which 
displays weak ferromagnetism (wf). The wf behaviour 
indicates that the ordered moments acquire a tiny ferro- 
magnetic component by a tiny tilt out of the layer plane. 
LaMnOa is also the parent of a family of magnetic mate- 
rials, the manganites, which are of current wide interest 
because of the interconnection between their structural, 
transport and magnetic properties. The features of this 
family, including colossal magnetoresistance (CMR) and 
a very large spin polarization of the conduction band 
around 1/3 doping, make them promising for future ap- 
plications. Moreover they are of prominent interest for 
the understanding of correlated electron systems. 

Manganites are characterized by the coexistence of sev- 
eral interactions on comparable energy scales, which pro- 
duce a complex interplay of magnetism, electron-lattice 
coupling, orbital and charge ordering. This shows up as 
soon as the system starts to be doped, typically by het- 
erovalent cation substitution, giving rise to phase separa- 
tion, magnetic polarons and eventually to the CMR phe- 
nomenon. However many of the basic ingredients, such as 
for instance orbital ordering and a large Jahn- Teller ( JT) 
distortion, are already present in pure LaMn03. Some 
authors actually question its description as a simple in- 
sulating antifcrromagnet, and suggest that the unusual 
layered magnetic structure is an effect of the static co- 
operative JT deformations on the magnetic interactions, 
as a sort of precursor to the double-exchange ferromag- 
netism of the doped compound 5 . These features under- 
lines the relevance of an accurate determination of the 



properties of the end member. 

Our single crystal investigation by Muon Spin Rotation 
(fiSlZ) attains two main results: we obtain a direct mi- 
croscopic determination of a few magnetic properties of 
this important material and we identify with a high con- 
fidence level the interstitial sites occupied by the muons. 
The first aspect covers static and dynamic critical expo- 
nents of the antiferromagnetic transition, as well as the 
initial process of weak ferromagnetic saturation. Critical 
phenomena have been addressed also in the substituted 
compounds, both by neutron scattering 1,2 and by macro- 
scopic magnetization^, hence it is of primary interest to 
obtain accurate measurements of the same quantities and 
behaviors with a microscopic probe in the end member, 
much less prone to stoichiomctry problem. Furthermore 
non trivial consequences of the weak ferromagnetism are 
also revealed. The second aspect is also of general rele- 
vance, since it contributes to the understanding of muon 
localization in crystals, with something more than a new 
single case. The insight by fiSlZ is often greatly enhanced 
by the knowledge of the muon site and the systematics 
in this field is far from complete. We believe that af- 
ter our contribution the situation in LaMnOa is better 
understood than, for instance, in all the cuprates. Hope- 
fully our results can be extended in part to other oxide 
and fluorite perovskites, including KCUF3 and the high 
T c materials. 

The paper is organized as follows: Sec. ITT1 briefly re- 
views the structural and magnetic properties of the com- 
pound, some details about our sample and a brief techni- 
cal introduction to [iSlZ; Sec. IIIII illustrates our results, 
including the identification of the muon stopping sites 
based on dipolar sum calculations and on a refined fitting 
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FIG. 1: Orthorhombic cell of LaMn03 with the A-type mag- 
netic structure (distortions and the wf canting are not shown 
for clarity). 



procedure, the zero field determination of the staggered 
magnetization, the dynamical critical behavior and the 
saturation of the weak ferromagnetic domain structure; 
These results are discussed in sec. IIVI Finally details on 
the dipolar sums and on the refined global fit are given 
in the appendices. 



II. EXPERIMENTAL 

A. LaMn0 3 

The crystal structure of LaMnC>3 is orthorhombic, be- 
longing to space group Pnmb, with lattice parameters^ 
a = 5.7391(2) A, b = 5.5319(2) A, c = 7.6721(2) A at 
14 K. The unit orthorhombic cell is shown in fig. ^ 
together with the pseudocubic cell (dashed). In the fol- 
lowing the subscripts 0>p refer to these two cells, respec- 
tively. Also shown is the magnetic structure determined 
by neutron diffraction below T/v = 139.5 K. The spin is 
along s = [010] = [110] p and the periodicity, given by 
qAF = (00i)p = (001) o , produces a staggered stacking 
of ferromagnetic planes (referred to as A-type antiferro- 
magnetism) . 

The absence of the inversion symmetry around the 
midpoint between neighboring Mn ions along [001] im- 
plies the presence of an antisymmetric exchange, the so 
called Dzialoshinski-Moriya interaction^^, which pro- 
duces weak ferromagnetism, i.e. a tiny canting of the 
spin from the orientation of a pure collincar antifcr- 
romagnet. This is indeed experimentally determined 
by magnetization^ and antiferromagnetic resonance^, 
with a canting angle of 2 degrees, while neutron diffrac- 
tion datai 3 - set an upper limit of 1 degree to its value. 
This second determination is insensitive to any accidental 
small fraction of the ferromagnetic compositions, hence 
it should reflect the intrinsic local value. 

Our sample is a microtwinned single crystal grown by 
the float zone method^ at the Laboratoire de Chimic 
Physique, Universite Paris-Sud, cut into a slice of roughly 
24x7x2 mm 3 . Six families of cubic microtwins are 



present, obtained one from another by interchange of 
±[100]p,±[010] p and ±[001] p . 



B. The /J.S7Z experiment 

The experiments were performed on the GPS 
instrument 15 of the Paul Scherrer Institut, Villigen 
(Switzerland). The normal to the largest surface of the 
crystal - close to a [111] direction - was oriented par- 
allel to the muon beam direction, along the laboratory 
frame x axis. The veto counter setup (VCS) was used to 
minimize positron counts from muons missing the small 
samplei 5 -. 

In a fiSlZ experiment the positron counts from im- 
planted ^ + decays reflect the spatial asymmetry of the 
decay probability, peaked around the initial \i spin di- 
rection ( — J(0)/|J|. Hence the count rate N(t) = 
N exp(— t/r)[l + A G(t)] is modulated by the muon 
polarization, i.e. the self correlation function of the 
muon spin projection in the counter direction £, G(t) = 
G^(t) — (I^(0)I^(t)), averaged over the implanted en- 
semble. Here r = 2.197/xs is the muon lifetime, Nq the 
initial positron count rate and ^4^ (£) = AoG(t) is often 
referred to as the muon experimental asymmetry func- 
tion. 

For a coherent precession around a transverse magnetic 
field B ± C, (e.g. in a calibration T > Tm experiment) 
one has G(t) = cos(27T7 \B\ t) and one can thus measure 
the asymmetry parameter Aq (typically around 0.26 for 
each GPS detector). A damped harmonic precession is 
obtained also in magnetically ordered materials, where a 
local field B at the muon site is produced by the electron 
magnetic moments. It has a dipolar contribution, Bdi P , 
and a hyperfine contribution, Bhf- The former is: 

n j n 

where x][ 3 with «,j = x,y,z are the Cartesian coordi- 
nates and r n the modulus of the vector joining the n-th 
magnetic lattice site (with magnetic moment Qf-iB^n 

) to 

the muon site, r n = R n — r^. The hyperfine field may be 
written as Bhf — Y^ j C J ■ Sj with Sj the spin of a near- 
est neighbor magnetic ion, the hyperfine path in oxides 
being through oxygen ions. 

Muons may occupy inequivalent sites. In the ordered 
phase and in zero external field the presence of cubic 
microtwins reproduces the usual powder average: 

Mt) = E^/ t cos^[2cos(2^ 7 |B^)e- t/T2 - 

i 

+ e"*/ T »] (2) 

where is the fraction of the muon ensemble localized 
at a site characterized by the spontaneous local field _B; ; 
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0£ = cos -1 (C - is the angle between the initial muon spin 
direction and the detector axis; the transverse (precess- 
ing) component is subject to spin-spin relaxation with 
rate T 2 _1 , while the longitudinal component is subject to 
spin-lattice relaxation with rate T± . 

In GPS the positron detectors lie in the y plane (with 
£ = ±z or ±i) and the cryostat lies along the y direction. 
A small magnetic field (5 mT) can be applied along y for 
calibration; larger fields may be applied along x. The 
initial muon spin direction lies in the y plane at a chosen 
angle 170 < ip — cos _1 (C • x) < 130 deg, depending on 
the settings of a cross-field device (spin rotator) in the 
beamline. The spin rotator was generally on (ip = 130 
deg) to optimize the efficiency in the VCS. 

Data analysis was performed by means of the MuZen 
suited. The basic approach to data analysis was by 
Fourier transforming or fitting to specific models an aver- 
aged muon asymmetry function, experimentally obtained 
from opposite groups of detectors. E.g. with ip — 130 deg 
(spin rotator on) 



A(t) 



N x + N z -a(N 1F + N-) 
N x + N z + a(Nx+N-) 



(3) 



where a is an average correction factor for the different 
detector efficiencies. This approach loses some of the 
phase information and the average over a large solid angle 
reduces the observed muon asymmetry, but it allows the 
implementation of the complete set of fitting functions for 
A(t) available in MuZen. A more powerful simultaneous 
fit of the four detectors, described in sec. El was used for 
fit refinements. 



III. RESULTS 

Muon polarization data at three temperatures below 
Tjv are shown in fig. together with their best fit (solid 
curves). Two precession frequencies are found at each 
temperature, indicating two magnetically inequivalent 
stopping sites for the muon. 



A. The muon sites 

The site identification is derived by comparing an anal- 
ysis based on dipolar sums with the results of two ex- 
perimental features of the internal fields probed by the 
muon: their magnitude, extrapolated to zero tempera- 
ture - -B° = 0.631(1) T and B% = 0.954(1) T - and their 
direction. The magnitude is measured directly by the 
precession frequency, while the direction is revealed by 
the frequency splitting upon application of an external 
field. Since the sample is twinned there was no scope 
for a full scan of the angular dependence of the resulting 
frequencies. However the most important symmetry is 
already apparent from fig. [3J where the two top curves 
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FIG. 2: Muon polarization, as from eq. at three tem- 
peratures below TV, displaced vertically for clarity; the solid 
curves are the best fits according to eq. El 



compare the Fourier transform (FT) of the muon po- 
larization signal with and without an external field ap- 
plied close to the [lll] p direction. In fig. |3 spectrum 
a) shows a low frequency precession around a local field 
B x = 0.5339(1) T at 70 K, which shifts to higher fre- 
quency under an external field of 100 mT (spectrum b). 
At the same time the high frequency precession around 
B 2 = 0.8150(2) T in spectrum a) undergoes a splitting in 
spectrum b). The inset shows a geometry which justifies 
this behaviour: the lower frequency corresponds to an in- 
ternal field B\ along [001] (the wf axis) while the higher 
frequencies correspond to an internal field perpendicular 
to the wf axis. The external field saturates the sample 
and adds vectorially as shown, giving rise to one lower 
and two higher precession frequencies. The geometry and 
magnitude of the internal fields is discussed below. 



1. Low field muon site 

We start by noting that an identical behavior is de- 
tected in several (iSIZ experiments performed on ortho- 
ferrite single crystals^. There the upward shift of the low 
frequency precession is justified on the basis of the fol- 
lowing two assumptions: a) a local field entirely of dipo- 
lar nature, and b) a muon site within the mirror plane 
z = 1/4 (or 3/4) of the orthorhombic cell. These assump- 
tions turn out to be very reasonable for the orthoferrites, 
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FIG. 3: Fourier transform of polarization, simulated from the 
best fit parameters for different crystal orientations, identi- 
fied by Euler angles (6, cj>, Q of the pseudocubic axes in the 
GPS reference frame. Plots are shifted both vertically and 
horizontally for clarity. 
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FIG. 4: Fourier transform of the best fit polarization at H = 
f00 mT, for different crystal orientations identified by Euler 
angles (6, of the pseudocubic axes in the GPS reference 
frame. Plots are shifted both vertically and horizontally for 
clarity. 



which are isostructural to LaMnOs, therefore we start 
from the same conjecture. Furthermore the assumption 
of Holzschuh et alH is justified by remarking that actu- 
ally point a) is a consequence of point b). Since the muon 
hyperfine coupling C is transferred through oxygen and 
also the relevant oxygen ions belong to the mirror plane, 
Bhf vanishes by symmetry (the hyperfine contributions 
from two Mn ions belonging to opposite sublattices can- 
cel). 

Appendix 1X1 discusses briefly the assignment: it shows 
that the non-vanishing dipolar field Bdip is along [001] p 
and it explains why a sufficiently large external field pro- 
duces only the upward frequency shift for such a site, 
which is a consequence of the weak ferromagnetism of 
LaMn0 3 . 

Notice that in a saturated sample, for a generic orien- 
tation of the twinned crystal, the local field at this kind 
of site 

= B dlp + fi H (4) 

should give rise to three inequivalent vector compositions 
(one per pair of twins sharing their c axis direction) , while 
only one upward shifted frequency is experimentally ob- 
served. The further reduction of multiplicity is due to the 
chosen crystal orientation which brings the three peaks 
in near coincidence. 

By contrast fig. 0] shows several simulations for more 
generic crystal orientations, where many frequency peaks 
are evident. Furthermore a fine structure is already ap- 
parent in all the peaks of fig. [21 curve b; this structure 



indicates that B is not exactly parallel to [lll] p in this 
experiment and the misalignment is actually resolved by 
the global best fit described in sec. [51 

In order to corroborate our assumptions for the low 
field site we performed a numerical evaluation of the dipo- 
lar sums (eq. ^| , considering the low temperature crystal 
structure^! and a moment of 4/x# per Mn ior£. We eval- 
uated the sums over a large spherical domain centered 
at the muon site, with a numerical error of less than 0.5 
mT. 

Fig. shows the magnitude of \Bdip(x, y)\ as a func- 
tion of the muon position (x,y) in the (001) o plane at 
z = 1/4; the large solid circles correspond to a radius of 
1.1 A centered on the O ions, equal to the average /^-0 
distance for oxidesi^, hence the muon site is expected to 
lie on its boundary. Candidate sites lie at the intersection 
of these curves with the field surface at the experimental 
value of 0.631(1) T. The inset shows their position in the 
mirror plane. 

One site, [0.389 0.937 0.25], is marked by a darker 
symbol, together with its symmetry replica, [0.889 0.563 
0.25]. It very closely corresponds to the main site iden- 
tified in the orthoferrites (see fig. 10, site 2, in refiii): 
they are farthest away from positive La ions, facing a 
large empty space and their cell coordinates also roughly 
agree. The correspondence implies that the internal field 
at this muon site scales with both lattice parameters and 
transition metal magnetic moments as dictated by the 
dipolar coupling. Hence the combined orthoferrite and 
LaMnOa experiments corroborate that the hyperfine con- 




FIG. 5: as a function of the muon position in [001] p 

mirror plane at z = 1/4; The large shaded circles of radius 
1.1A are centered at O sites; The inset shows the mirror plane 
with O (open circles) , La (filled circles) and candidate mu sites 
(open symbols, H sites are filled). 

tributions must be negligible. We shall refer to these and 
their replica in the z = 3/4 mirror plane as the Holzschuh 
(H) site. 

Holzschuh et al. calculated that also along the [i z) p 
line Bdip has the wanted field orientation. However along 
this locus a consistent hyperfine contribution is required 
to agree with the experimental internal field value. Since 
the hyperfine term in general is not collinear to the dipo- 
lar term, the required experimental symmetry would be 
lost. The TL site alone can satisfy both magnitude and 
direction constraints. 

The exact vector composition of HqH and Bdi P for 
the experiment of fig. OH were obtained from a global fit 
of the [iSlZ data - described in detail in appendix [B] 
which yields the following /iSTZ assessment of the crys- 
tal orientation: the three pseudocubic axes had director 
cosines of [0.74(3) 0.51(3) 0.43(4)] with the external field. 
This is in reasonable agreement with a subsequent neu- 
tron diffraction determination [0.75(6) 0.57(6) 0.36(6)], 
where the large error bars come mainly from systematic 
errors of ~ 4 deg in the sample realignment between the 
two experiments. 

2. High field muon site 

Fig. 13b reveals that the application of the external 
field in the chosen experimental geometry gives rise to 
two precessions equally separated on either sides of the 
zero field high frequency peak. Thus i) the high fre- 
quency site does not belong to the mirror plane and ii) 
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FIG. 6: The basal plane a) in the undistorted cell; b) when an 
O-Mu-0 center is formed in C. Mn ions (filled circles); O ions 
(open circles) are centered 0.03 A above - solid - or below - 
dashed - the plane; the C site (dotted circle); the points on 
the solid squares are also loci where Bdi P + Bhf is parallel to 
[100] o (see text). 



only two main peaks are observed. This implies again a 
highly symmetric site since the generic site in a Pmnb 
cell (type d in Wyckoff notation^, with multiplicity 8) 
in a generic orientation of the twinned crystal in the ex- 
ternal field would give rise, according to eq. Qi, to 48 
inequivalent vector compositions! Their number reduces 
greatly thanks to the chosen crystal orientation, but the 
experimental results are compatible only with the net in- 
ternal field B ( ii P +Bhf along one of the pseudocubic axes. 
This is shown in fig. 0] where Fourier spectra of the best 
fit function are simulated for different orientations of the 
crystal in the external field. It is clear that both the high 
symmetry of the muon site and the best-fit crystal ori- 
entation are rather stringent conditions: if the two were 
not simultaneously met the spectrum would be spread 
out into many subpeaks, mostly below noise level. We 
therefore conclude that the internal field Bdi P + Bhf lies 
along the local [110] p =[100] o direction, i.e. the direction 
of the ordered moments. 

Assuming isotropic hf coupling, both Bdip and Bhf 
must be parallel to [100] . This particular field orienta- 
tion is obtained along the following lines in the Pnmb 
cell: [x 0] o , [0 y 0] o , [x 0.5 0] o , [0.5 y Q] OJ [0.5 z] , and 
[0.5 0.5 z] , plus those equivalent by symmetry. Most of 
them coincide with the oxygen-cation bonds, thus being 
unlikely locations for a positive muon. However the [x 

0. 5 0] o and [0.5 y 0] o lines connect two magnetic ions and 
intersect at the center of the basal plane, in C=[0.5 0.5 

01 . The region around C is empty and it is a possible lo- 
cation for the high frequency muon site. Site C is shown 
in fig. [5p,, together with a few of the mentioned lines. 

There is an additional reason to consider site C in par- 
ticular. It is well known that muons implanted in several 
fluorides^ end up in a F-Mu-F structure, where two F 
ions are drawn towards the muon sitting at an interstice 
midway between them, in a muon version of the hydrogen 
bond. The bond distance in this center is remarkably the 
same (1.16 A) throughout a large series of compounds. 
We have recently observed^ the same center in KCUF3, 
which is essentially isostructural to LaMnC>3. We sug- 
gest that an analogous O-Mu-0 center could be formed 
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FIG. 7: Temperature dependence of the Ti and C internal 
fields. The solid lines are the best fit of the critical behavior. 
The inset replots the data on a log-log scale; the solid lines 
are the same best fits. 



here, exactly at site C, as shown in fig. [SJd. The resulting 
overlap with O could easily justify the required hyperfine 
contribution of Bhf = 0.51 T, in addition to the calcu- 
lated B dlp {C) = 0.41 T. 

Finally one may note from the inset of fig. [S]that the 
O ions in the mirror planes cannot be drawn together be- 
cause of the La ion hindrance, whence the Ti site remains 
locally favored in that plane. 



B. Critical behaviour 

The zero field experiments provide a measurement of 
the temperature dependence of the staggered magnetiza- 
tion. Figure [7| shows the behavior of the internal fields 
at both muon sites. They are proportional to the local 
magnetic moment on Mn, i.e. to the magnetic order pa- 
rameter, by means of the coupling to the muon spin, as 
shown e.g. in eq. ^ Hence we may extract the critical 
exponent of the order parameter, /3, from them; to this 
end they are replotted as a function of reduced temper- 
ature £ = 1 — T/Tn on a log-log scale in the inset of fig. 
|7| The data follow a power law e° on a remarkably large 
temperature interval, at least up to e — 0.6. The crit- 
ical exponent extracted from this interval, (3 = 0.36(2) 
is that predicted by the 3D Heisenberg model^, within 
experimental errors. 

Incidentally, this analysis reveals a perfect agreement 
between the two muon sites, which confirms that their 
occupancy is stable in the observed temperature interval. 

We now turn to the muon relaxations which give access 
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FIG. 8: Critical divergence of the transverse and longitudinal 
relaxation rates; circles: T 2 -1 (open - Ti site, filled - C site), 
solid triangle T^ 1 . The inset shows the critical part on a 
log- log plot. 



to the dynamical aspects of the critical behaviour 24 ' 25 26 . 
Approaching T/y from below fluctuations the AF order 
appear. The closer one gets to TV, the larger are the 
fluctuating clusters and the slower is the rate of spin 
fluctuation, both effects contributing to a diverging spin 
relaxation. 

In preliminary reports2i2& we wrongly remarked the 
absence of a critical divergence of the relaxation rates in 
LaMnC>3. We had overlooked that in those early data 
the sample had experienced an external field at low tem- 
perature, resulting in inhomogeneous broadening effects 
which masked the critical divergence. Figure [S] shows 
instead that in proper zero field experiments both T^ 1 
and T 2 _1 diverge at Tjv as e~ n . By subtracting the non- 
critical contribution (indicated by the dashed lines in the 
figure) a dynamical critical exponent n = 0.7(1) may be 
extracted from the relaxation rates^ 4 *^ 5 ^ 2 ^, which is in 
agreement with the theoretical best estimates^ based on 
the Ising model (n sa 2/3 - cfr. with Heisenberg model, 
n Ri 0.35). This is indeed what is expected, since the 
anisotropy of the interaction dominates in the critical 
regime^ .. 



C. Weak ferromagnetic domain structure 

Fig. shows the FT spectra of several data sets 
recorded in increasing external longitudinal fields (1 < 
Mo-ff < 100 mT, apodization exponential filter of 100 
kHz) and figEK displays the field dependence of the cor- 
responding frequencies. 

Above 20 mT the features described in details in sec. 
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FIG. 9: Longitudinal field data: a) Frequencies obtained 
from the global fit (crosses) and from the FT (open circles); 
the lines are guides to highlight the low (dotted) and high 
(dashed) field behaviour; b) FT spectra of the data (displaced 
diagonally for clarity). 

IIII Al are apparent from both plots: the frequency belong- 
ing to the Ti. site shifts upwards and eventually gives rise 
to three satellites, while the C site gives rise to two groups 
of satellites on either sides of its zero field resonance. 

The extrapolation of this high field behaviour to low 
fields evidences an effective saturation induction, /iqHs = 
7 ± 1 mT. This is best seen from the average field de- 
pendence of the 7i peaks (although it holds also for the 
C peaks). The dashed line is optimized on the high field 
values of the central fit frequency, while the low field data 
agree better with the constant dotted line. The quoted 
demagnetization is taken from the intercept of the two 
lines. This behaviour resembles that of a soft magnet, 
although LaMnOs is definitely not soft: the total mag- 
netic energy, oc B ■ H is minimized for B — 0, hence 
external fields Hq < H s are shielded inside the sample. 



IV. DISCUSSION AND CONCLUSIONS 

We have definitively confirmed the TC site location for 
pseudocubic perovskites first suggested by Holzschuh et 
alw^ in orthofcrrites. This indicated that the analogy 
between manganites and orthoferrites is worth pursuing. 
Holzschuh et al. detected three magnetically inequivalent 
sites for the muon. The missing fraction in our experi- 



ment probably justifies a third site, with either a very 
large internal field (a conservative estimate of the detec- 
tion frequency-passband cut-off is 400 MHz) or very large 
relaxation rates. Holzschuh et al. discuss the metastabil- 
ity of their high frequency sites, which depopulate above 
w 150 K in favor of the low frequency H site. We cannot 
directly confirm this feature since Tjy in LaMnOs is too 
low and the two observed sites are both stable up to that 
temperature. 

We have found a strong analogy between the high field 
site and the F-Mu-F center detected in fluorides^. As a 
matter of fact the formation of an O-Mu-0 center, the 
muon version of the hydrogen-bonding, is a very natural 
expectation for partially ionic crystals. 

After the first draft of the present paper we became 
aware of a preprint^ on [iSlZ in Lai-^Ca^MnOs pow- 
ders with < x < 1. They confirm our earlier reports 
of two sites2L2& in LaMnC>3 and agree on the TL site at- 
tribution. Their data, however, are consistent with two 
sites only, while they observe three sites in the other end- 
member, CaMnC>3. 

Let us now discuss briefly the critical behaviour as seen 
from the local fields and from the relaxations. The log- 
log plots of fig. [7| (inset) shows a remarkably large range 
on which the behaviour holds, and the best fit for 
both sets yields with T N = 139.3(1) K and (3 = 0.36(2). 
A cross-over from Heisenberg (/? = 0.365(3)) to 3D-Ising 
behaviour— ((3 = 0.325(2)), is expected close enough to 
Tn, since the anisotropy which is present in the ordered 
phase must eventually dominate, for very small e. This 
cross-over is not directly observed in our magnetization 
data, but this is hardly significative, in view of the small 
variation of (3 between the two limits. The Ising-type 
behaviour does indeed show up in the relaxation rates 
of fig. 03 whose critical exponent (related— to the dy- 
namical exponent z) displays a value, n — 0.7(1), in full 
agreement with 3D-Ising predictions. 

We must mention that the muon determination of the 
temperature dependence of the staggered magnetization, 
M M (T), does not fully agree with that obtained by neu- 
tron diffraction, M n (T), on the same crystal^. Assum- 
ing that the two data sets agree at T = 0, where the 
muon site attribution is fully consistent with w 4|ig on 
Mn determined by neutrons, the deviation - M M < M n 
- increases with temperature. Although the origin of 
this discrepancy is not understood, we assume that both 
techniques measure intrinsic sample properties. We hope 
that our planned Mn NMR experiments will help clari- 
fying this issue. 

Let us finally discuss the observed wf saturation. For 
an ideal soft wire one would have H s = M s . Since the 
sample is twinned we must expect a shielding field from 
the macroscopic average over a large number of cubic 
twins to mimic that of a soft material, albeit with sat- 
uration magnetization M s /3. This reduction factor is 
most clearly explained thinking of a twinned single crys- 
tal with Ho || [001], the easy wf axis for one family 
of twins: only their contribution, 1/3 of the total, pro- 
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duce a shielding. Since LaMnC>3 has a large anisotropy, 
only the projection of the external field H along the 
easy axes may be shielded by demagnetization, hence 
H s cos9 = M s /3, where 9 is the angle between the 
easy wf axis and H (in our experiment the average 
of cos6> is 0.53). Approximating the sample shape to 
an ellipsoid, with demagnetizing factor N, one finally 
has that H s = NM S / (3 cos 9). An estimated value of 
N = 0.17 may be derived from standard formulas^!, to 
yield M s = 3.4(6) kA/m, i.e. 0.22(4) fi B at 90 K, consis- 
tent with magnetization measurement*^. This implies a 
canting angle of order 2 degrees, larger than the 1 deg up- 
per limit set by neutrons^. However our determination 
is essentially equivalent to a magnetization measurement 
and it is not a true microscopic determination, although 
it is detected via a local probe. 

A closer inspection at the 10 and 20 mT spectra of 
fig. El is in qualitative agreement with this picture, with 
a small further coercive effect: at 10 mT an incipient 
splitting of the TL frequency is apparent, with a negligible 
domain polarization, while at 20 mT a low frequency 
shoulder of the Ti. peak still appears, from the minority 
domains. 

In conclusion we have clarified the identification of the 
two muon sites in LaMnOa, the TC and C sites. Our single 
crystal data provide a definitive assessment of the former, 
while the latter is proposed here for the first time. The 
C site, analogous to what is observed in many fluorides, 
might be more frequent than previously thought. These 
assignments allow us to extract an estimate of the wf 
canting angle. The critical behaviour of this material 
agrees with its current view as a Heisenberg system with 
a cross-over to Ising behaviour due to anisotropy. 



APPENDIX A: DIPOLAR SUMS 

Let us consider the dipolar contribution to local field 
at a muon site, given by eq. ^ ln the ordered state one 
may replace all with either S 3 a or Sp, depending on 
whether n belongs to sublattice a or /3, to yield 



Mn 



B!, 



u a p 



(Al) 



where the two tensors D represent the sum of eq. Q 
restricted to each sublattice. For a collinear antiferro- 
magnet (a very good approximation for LaMnOa in view 
of the minute canting) Sp — —S a = S, hence B = 
(D a — Dp) ■ S and one may refer simply to D — D a — Dp. 

Let us consider first the low frequency site. Figure 
1101 graphically shows that the total dipolar tensor D for 
muon sites in the mirror plane z — 1/4, 3/4 of the Pnmb 
cell possesses the following symmetry: 
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d 3 
d 2 
d 3 d 2 
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FIG. 10: Symmetries of the local dipolar tensor at equiva- 
lent muon sites in the mirror planes z — 1/4, 3/4; three Mn 
moments along c are shown (open and grey circles with thick 
spin arrows) with their dipole field lines (dashed) through the 
muon sites. Other arrows represent the local magnetic field 
at the muon from a pair of symmetric moments: a) Mn spin 
component along c; b) Mn spin component normal to c, in the 
plane containing the muon sites; The third component (not 
shown) gives vanishing contribution by symmetry. 



which guarantees that, for s = [110] p , Bdip lies along 
[001] p . 

If we neglect the tiny spin canting the Mn moments lie 
in the [001] plane, hence we can concentrate on fig. HObh 
which illustrates the relevant case for both mirror planes. 
Since the two Mn sublattices are symmetric with respect 
to the mirror plane, the dipolar field properties are iden- 
tified considering just two equivalent Mn ions, one per 
sublattice. Only the components of the Mn magnetic mo- 
ments in the plane of the drawing are considered, since 
the perpendicular components give opposite dipolar field 
contributions which cancel. The thinner arrows in the 
figure indicate the dipolar term from each Mn ion at a 
generic site in the mirror plane: it is graphically evident 
that their sum lies along [001], as predicted by eq. IA2I 
Similar arguments hold for the other components of the 
tensor. 

Notice that the figure shows just one AF domain (the 
bottom Mn moment belongs, say, to the a sublattice); 
the opposite direction of Bdi P is obtained in the other 
AF domain (where the bottom Mn moment belongs to 
the (3 sublattice; imagine the figure with all arrows re- 
versed). However according to the Pnmb symmetry the 
[3 domain must also have its wf component in the oppo- 
site direction along [001] p . 

Based upon these consideration we may understand 
the field behavior of the fiSTZ precessions. The wf spin 
component is totally negligible in our dipolar sums, but 
it gives a net wf moment M w * to each domain, which 
couples to an external field H. When H saturates the 
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macroscopic wf moment only one of the two domains 
survives - specifically that with H ■ M w ^ > 0. It ap- 
pears that M w * and Bdi P at the 7i muon site are point- 
ing in the same direction, so that when the component 
of H along [001] is larger than the saturation value one 
has H ■ Bdip > 0, hence the precession frequency shifts 
upwards, as in fig. [3] By contrast at the C site the in- 
ternal field is in the direction of the spin and its the vec- 
tor composition with /j,qH is not influenced by the AF 
domain structure (all possible geometries are obtained 
within each domain). 



APPENDIX B: GLOBAL FIT 

We have employed a procedure to fit the whole set of 
data recorded at each temperature and field. The global 
chi square, X 2 = ~ N° h (t i )) 2 /N a (t i ), is ob- 

tained from the count rates N a of the four detectors di- 
rected along a = ±x, ±y, at each recorded time £j. The 
model function Nth, briefly described in sec. Ill Bl de- 
pends on several parameters, of which only few - field 
intensities, their orientations and the relaxation rates - 
are informative of the magnetic properties of LaMn03. 
The rest are of lesser interest in this respect, such as the 
angle ip — cos _1 (C • z) defining the initial muon spin di- 
rection, the experimental asymmetry A of the reference x 
detector, the muon fractions, fi and the individual initial 
count rates, Nff, instrumental delay times, t a and rela- 
tive efficiencies e a for all the a detectors. The asymmetry 
parameter of each detector is Afce 01 , and the efficiency 
of the reference detector is fixed, e x = 1. All these are 
a byproduct of the fit procedure, and a number of them 
may be accurately determined in separate high statistics 
calibration experiments. 

All models generally include one or more contribution 
to the polarization functions, as in the case described by 
eq. |2 which points out that quite generally each local 
field Bi identifies a muon fraction fi characterized by its 
transverse and longitudinal relaxation rates, \jTn and 
1/Tu. When fitting more than one fraction, however, it 
is often impossible to distinguish individual longitudinal 
rates and a single \jT\ value is employed. 

The advantage of the global approach is that the direc- 
tion of the local field Bi may be extracted from the data, 
since it determines the aperture of the muon spin preces- 
sion cone. This is incidentally not true for the zero-field 
sum over our cubic twins, discussed in eq. |2 since the 
information on the local field direction is lost by sym- 
metry. Nonetheless this zero field model is very easy to 
implement in the global fit scheme and the field values of 
section All BI are indeed obtained from such a fit. 

The global approach becomes more valuable in an ap- 
plied external field H , were the direction of the internal 
field Bhf may be determined. In this case the direction of 
the local field B^f + hqH in each twin may be expressed 
in terms of the three Euler angles (0, <fi, £) which define 
the crystal orientation in the laboratory xyz frame. 
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FIG. 11: Some parameters of the zero field global fit: muon 
fractions at the two sites and initial muon phase. The inset 
compares the detector time delays in zero field (open symbols) 
and in longitudinal fields (filled symbols). 



Therefore, following the discussion of appendix lAl we 
have implemented two further models: one to account for 
the 7i site in the saturated wf state (spontaneous field 
B,iip along the local c direction, in three twin replicae, 
chosen with the condition that Bdi P -H > 0); another one 
to account for the C site (spontaneous field Bhf along the 
local electron moment direction - [010] o - in six distin- 
guishable replicae). The analytic expression of the corre- 
sponding functions N? h (t) , involving their dependence on 
(8, <fr, and the vector composition of the internal and 
external fields, was generated and checked by means of a 
simple Maple V program. 

The best-global-fit values discussed in section IlII Al re- 
fer to a model function including both sites. Because 
of the relatively large number of parameters the conver- 
gence of the global fit requires some practice, but reliable 
fits were always obtained with y 2 values better than 1.1 
per degree of freedom. 

The fit quality is good but not perfect, as it may be 
judged from the comparison between the FT of the data 
and of the fitting function in fig. [3]d,c. The small dis- 
agreement may be imputed to simplifying assumptions in 
the fit function, in which the crystal twins are weighted 
equally, while x-ray diffraction in similar samples showed 
slight preponderance of one twinAi, and relaxations are 
all assumed Lorentzian, while probably a Gaussian Ti 
rate is more appropriate. 

The behaviour of some of the less interesting fit pa- 
rameters confirms the soundness of the model. Figure ITT1 
shows the temperature dependence of the muon fractions 
and of the initial muon spin phase ip, obtained from the 
zero field data. The phase ip is perfectly stable, while the 
two fractions do show some small spurious correlation ef- 
fects close to the transition, where the relaxation rates 
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diverge. This is acceptable since it has no appreciable 
influence on the critical analysis of sec. 1111 Ul The inset 
shows the scattering of the time delays of the four de- 
tectors both in the zero field data (open symbols) and in 



the longitudinal field data, which are fitted with different 
functions. These parameters are stable and insensitive to 
both field and model function employed. 
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